Bose-Einstein condensation and superfluidity of magnetobiexcitons in quantum wefls' 

and graphene superlattices 
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We propose the Bose-Einstein condensation (BEC) and superfluidity of quasi-two-dimensional 
(2D) spatially indirect magnetobiexcitons in a slab of superlattice with alternating electron and 
hole layers consisting from the semiconducting quantum wells (QWs) and graphene superlattice 
in high magnetic field. The two different Hamiltonians of a dilute gas of magnetoexcitons with a 
dipole-dipole repulsion in superlattices consisting of both QWs and graphene layers in the limit of 
high magnetic field have been reduced to one effective Hamiltonian a dilute gas of two-dimensional 
excitons without magnetic field. Moreover, for A'^ excitons we have reduced the problem of 2 A*' x 2 
dimensional space onto the problem of A'^ x 2 dimensional space by integrating over the coordinates 
of the relative motion of an electron (e) and a hole (h). The instability of the ground state of 
the system of interacting two-dimensional indirect magnetoexcitons in a slab of superlattice with 
alternating electron and hole layers in high magnetic field is established. The stable system of indi- 
rect quasi-two-dimensional magnetobiexcitons, consisting of pair of indirect excitons with opposite 
dipole moments is considered. The density of superfiuid component ns{T) and the temperature of 
the Kosterlitz-Thouless phase transition to the superfiuid state in the system of two-dimensional 
indirect magnetobiexcitons, interacting as electrical quadrupoles, are obtained for both QW and 
graphene realizations. 

PACS numbers: 71.35.Ji, 71.35.Lk, 71.35.-y, 68.65.Cd, 73.21.Cd 



The many-particle systems of the spatially-indirect ex- 
citons in coupled quantum wells (CQWs) in the presence 
or absence of a magnetic field B have been the subject 
of recent experimental studies. ^'^''^ These systems are 
of interest, in particular, in connection with the possi- 
bility of BEC and superfluidity of indirect excitons or 
electron-hole pairs, which would manifest itself in the 
CQWs as persistent electrical currents in each well and 
also through coherent optical properties and Josephson 
phenomena^i^i^iii^i. In strong magnetic fields {B > 7T) 
2D excitons survive in a substantially wider temperature 
region, as the exciton binding energies increase with mag- 
netic field^ '-^°i-^-^'-^^i-^^'-^'*i^^ . The problem of essential inter- 
est is also collective properties of magnetoexcitons in high 
magnetic fields in superlattices and layered systemic. 

Recent technological advances have allowed the pro- 
duction of graphene, which is a 2D honeycomb lattice 
of carbon atoms that form the basic planar structure in 
graphiteiii^. Graphene has been attracting a great deal 
of experimental and theoretical attention because of un- 
usual properties in its bandstructureilii^iiSii^. It is a 
gapless semiconductor with massless electrons and holes 
which have been described as Dirac-fermions^. Since 
there is no gap between the conduction and valence bands 
in graphene without magnetic field, the screening effects 
result in the absence of excitons in graphene in the ab- 
sence of a magnetic field. A strong magnetic field pro- 
duces a gap since Landau levels form the discrete energy 
spectrum. The gap reduces screening and leads to the 
formation of magnetoexcitons. We consider magnetoex- 
citons in the superlattices with alternating electronic (e) 



and holes (h) QWs as well as graphene layers (GLs) and 
suppose that recombination times may be much greater 
than relaxation times Tr due to small overlapping of the 
spatially separation of e- and h- wave functions in QWs 
or GLs. In this case electrons and holes are characterized 
by different quasiequilibrium chemical potentials. There- 
fore, in the system of indirect excitons in superlattices 
the quasiequilibrium phases appear, as, for example, in 
CQWs^i^. While coupled- well structures with spatially 
separated electrons and holes are typically considered to 
be under applied electric field, which separates electrons 
and holes in different QWs — , we assume there are no 
external fields applied to a slab of superlattice. If elec- 
tron and hole QWs alternate, there are excitons with 
parallel dipole moments in one pair of wells, but dipole 
moments of excitons in other neighboring pairs of neigh- 
boring wells have opposite direction. This fact leads to 
essential distinction of properties oi e — h system in su- 
perlattices from one for CQWs with spatially separated 
electrons and holes, where indirect exciton system is sta- 
ble due to a dipole-dipole repulsion of all excitons. This 
difference manifests itself already beginning from three- 
layer e~h— eovh — e — h system. The superlattice of 
the GLs with the electrons and holes separated in alter- 
nating layers can be created by controlling the chemical 
potential of the charge carriers due to the doping of the 
GLs by the charged impurities^^ . 

In this Letter wc propose a new physical realization of 
magnetoexcitonic BEC and superfluidity in superlattices 
with alternating electronic and hole layers, that is in a 
sense representing an array of QWs or GLs with spatially 
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separated electrons and holes in high magnetic field. We 
reduce the problem of magnetoexcitons to the problem 
of excitons at i? = 0. The instability of the ground state 
of the system of interacting indirect excitons in slab of 
superlattice with alternating e— and h— layers is estab- 
lished in a strong magnetic field. Two-dimensional indi- 
rect magnetobiexcitons, consisting of the indirect mag- 
netoexcitons with opposite dipole moments, are consid- 
ered in a high magnetic field. These magnetobiexcitons 
repel as electrical quadrupoles at long distances. As 
a result, the system of indirect magnetobiexcitons be- 
comes stable. Below the radius and the binding energy 
of indirect magnetobiexciton are calculated. By apUy- 
ing the ladder approximation we consider a collective 
spectrum of the weakly interacting by the quadrupole 
law two-dimensional indirect magnetobiexcitons and cal- 
culate their superfluid density ns{T) in superlattices at 
low temperatures T . We analyze the dependence of the 
Kosterlitz-Thouless transition^'^ temperature and super- 
fluid density on magnetic field. 

The Hamiltonian of a single 2D magnetoexciton is 
given in Refs. for CQWs and in Ref. [H for GLs. 

A conserved quantity for an isolated electron-hole pair in 
magnetic field B is the exciton generalized momentum P 
defined as P = —ihVe — iKS/h + e/c{Af. — Ah) — e/c[B x 
(re — r/i)], for the Dirac equation in GLs^i as well as for 
the Schrodinger equation in CQWs^iii"^. Here re and 
Yh are electron and hole locations along QWs, respec- 
tively, D is the distance between electron and hole QWs, 
e is the charge of an electron, c is the speed of light 
and e is a dielectric constant. We find it is convenient 
to work with the symmetric gauge for a vector potential 
Ae(/t) ~ 1/2[B X re(/i)], where re(;j) is the location of an 
electron or a hole, respectively. 

The Hamiltonian of a single isolated magnetoexciton 
without any random field (T4(re) — Vh{Yh) — 0) is 
commutated with P, and hence they have the same 
the eigenf unctions, which have the following form (see 
Refs. 

*fep(R,r) = e{'^(P+ti^^^)+^^'^}$fc(P,r), (1) 

where 7 = (m/j — me)/(m/i -I- me) (wg and rrih are the 
masses of an electron and hole in QWs, and 7 = for 
GLs), $fc(P, r) is a function of internal coordinates r and 
the eigenvalue of the generalized momentum P, and k 
represents the quantum numbers of exciton internal mo- 
tion. The expression for $fe(P,r) for CQWs is given in 
Ref. [g'Ts'l, and in Ref. [24*]) for graphene layers. In high 
magnetic fields the magnetoexcitonic quantum numbers 
k = {n+,n_} for an electron in Landau level and a 
hole in level 

For l arge elec tron-hole separation D ^ rs, where 
rs — \/fic/ (eB), transitions between Landau levels due 
to the Coulomb electron-hole attraction can be neglected, 
if the following condition is valid for the magnetoex- 
citonic binding energy Ei,: Et, = e^/(ef,£') <C huJc = 
TieB(me + mh) / (2memhc) for QWs; Eb = 4e^/(eD) < 



hvp/rs for GLs, where vp = ^/3at/{2fi) is the Fermi 
velocity of electrons for a lattice constant a — 2. 566 A 
and the overlap integral between nearest carbon atoms 
t w 2.71eV^^. This corresponds to high magnetic field 
B, large interlayer separation D and large dielectric con- 
stant of the insulator layer between the GLs. At small 
densities nrg <^ 1 the system of indirect excitons at low 
temperatures is the two-dimensional dilute weakly non- 
ideal Bose gas with normal to wells dipole moments d 
in the ground state (d ~ eD, D is the interwell sepa- 
ration), increasing with the distance between wells D. 
In contrast to ordinary excitons, for a low-density spa- 
tially indirect magnetoexciton system the main contri- 
bution to the energy is originated from dipole-dipole in- 
teractions J7_ and C/+ of magnetoexcitons with opposite 
and parallel dipoles, respectively. Two parallel {+) and 
opposite (— ) dipoles in low-density system interact as 
[/_!_ = — J7_ = / eR^ , where R is the distance between 
dipoles along wells planes. We suppose that D/R <C 1 
and L/R 1 {L is the mean distance between dipoles 
normal to the wells). We consider the case, when the 
number of QWs or GLs k in superlattice is restricted 
k <^ {Dy/nn) ^, where n is exciton surface density, and 
this is valid for small k or for sufficiently low exciton 
density. 

Due to the orthonormality of the single magnetoexci- 
ton wave functions <I>„^,„„ (0, r) (r = re — r/i) the projec- 
tion of the the many-particle Hamiltonian onto the lowest 
Landau level results in the effective Hamiltonian, which 
does not reflect the spinor nature of the four-component 
magnetoexcitonic wave functions in graphene. Typically, 
the value of r is r^, and P h/vB- In this approxi- 
mation, the effective Hamiltonian Heg in the momentum 
representation P in the subspace the lowest Landau level 
(for QWs n+ = Ti- = 0; for GLs n+ = n_ = 1) has the 
same form (compare with Ref.Q) as for two-dimensional 
boson system without a magnetic field. Only differences 
are that, instead of the exciton mass M = me + rrih and 
ordinary momenta, we have the magnetoexciton mag- 
netic mass, which depends on B and D and magnetic mo- 
menta, respectively. For the lowest Landau level we de- 
note the spectrum of the single exciton eo(P) = eoo(P)- 
The Hamiltonian can be represented in the form: Htot = 
Hq + Hint , where Hq is the effective Hamiltonian of the 
system of noninteracting magnetoexcitons: 

^0-51 ^o(p)(a+ap + b+b^ + a+_^a_^ + b+_^b_^). (2) 
p 

In Eq. ^ £o{p) — /{^niB) is the spectrum of isolated 
two-dimensional indirect magnetoexciton; p represents 
the excitonic magnetic momentum and a+, a^, b^, b^ 
are creation and annihilation operators of magnetoexci- 
tons with up and down dipoles. For an isolated magne- 
toexciton on the lowest Landau level at the small mag- 
netic momenta under consideration, £o(P) ~ ^^/(Stib), 
where is the effective magnetic mass of a magnetoex- 
citon in the lowest Landau level and is a function of the 
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distance D between e - and h - layers and magnetic field 
B In a strong magnetic field at D ^ the exciton 
magnetic mass is mB{D) ~ eD^ /{e'^r%) for QWs"'^'^ and 
msiD) = eD3/(4e2r|) for GLs2i. The effective Hamil- 
tonian of the interaction between magnetoexcitons Hi„t 
is: 



Hint — 2S Spi+p: 



P3+P4 

n+ n+ h h 

P4 P3 P2 Pi 
^P4^P3^P2^Pl 



(flp, at^a 



P3"P2'3^Pl 
^P4^P3«"''" 



P2"'Pl 



(3) 



where S is the surface of the system. Let us consider the 
temperature T — Q. We apply Bogolubov approxima- 
tion and assume [N ~ Nq)/Nq ^ 1, where N and A^o are 
the total number of particles and the number of particles 
in condensate, respectively. In Bogolubov approxima- 
tion the interaction between non-condensate particles is 
neglected, and only the interactions between condensate 
particles and excited particles with condensate particles 
are considered. Therefore, in the total Hamiltonian the 
terms, arising from first and second terms of the Hamil- 
tonian ([3]), which describe the repulsion of the indirect 
magnetoexcitons with parallel dipole moments, are com- 
pensating by other terms of this Hamiltonian, describing 
the attraction of indirect magnetoexcitons with opposite 
dipoles. As the result only terms describing the attrac- 
tion survive. Let us diagonalize the total Hamiltonian by 
using of the Bogolubov type unitary transformation 



=("p 



Apl3+p - 



Bpa_p 



(4) 



where the coefficients Ap and Bp are found from the con- 
dition of vanishing of coefficients at nondiagonal terms 
in the Hamiltonian. As the result we obtain Htot — 
^p_^Qe(p)(Q!+ap + /3p/3p) with the spectrum of quasi- 

1 /2 

particles e{p): e{p) — [(eo(p))^ — (nU)'^] . At small 
momenta p < \/2mBnU the spectrum of excitations be- 
comes imaginary. Hence, the system of weakly inter- 
acting indirect magnetoexcitons in a slab of superlattice 
is instable. The instability of magnetoexcitons becomes 
stronger as magnetic field higher, because ms increases 
with the increment of magnetic field, and, therefore, the 
region of p resulting in the imaginary collective spectrum 
increases as B increases. 

This, on first view, strange result can be illustrated by 
the following example. There are equal number of dipoles 
oriented up and down. Let us consider four dipoles, two 
of them being oriented up and two — down. It is easy to 
count that number of repelling pairs is smaller than that 
of attracting ones. The prevailing of attraction leads to 
instability. 

We assume the energy degeneracy respect to two pos- 
sible spin projections in QWs and graphene and two 
graphene valleys (two pseudospins). Since electrons on 
a graphene lattice can be in two valleys, there are four 
types of excitons in bilayer graphene. Due to the fact 
that all these types of excitons have identical envelope 
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FIG. 1: 2D indirect magnetobiexcitons consisting of indirect 
magnetoexcitons with opposite dipole moments, located in 
neighboring pairs of QWs (GLs). 

wave functions and energies^i, we consider below only 
excitons in one valley. Also, we use no = n/{As) as the 
density of excitons in graphene superlattice, with n de- 
noting the total density of excitons, and s the spin degen- 
eracy, which equals to 4 for magnetoexcitons in bilayer 
graphene. Besides, we use tt-q = n/ s as the density of ex- 
citons in QWs, with the spin degeneracy s, which equals 
to 4 for magnetoexcitons in QWs. 

Let us consider as the ground state of the system the 
low-density weakly nonideal gas of two-dimensional indi- 
rect magnetobiexcitons, created by indirect magnetoex- 
citons with opposite dipoles in neighboring pairs of wells 
(Fig. [1]). The small parameter for the adiabatic approx- 
imation is the numerical small parameter which is equal 
to the ratio of magnetobiexciton and magnetoexciton en- 
ergies or the ratio between radii of magnetoexciton and 
magnetobiexciton along QWs or GLs. These parame- 
ters are small, and they are even smaller than analogous 
parameters for atoms and molecules. The smallness of 
these parameters will be verified below by the results 
of the calculation of indirect magnetobiexciton. Here it 
was assumed, that the distance between wells (GLs) D 
is greater than the radius of indirect magnetobiexciton 
flf). The potential energy of interaction between indirect 
magnetoexcitons with opposite dipoles U (r) has the form 
(r is the distance between indirect magnetoexcitons along 
QWs/or GLs): 



U{r) = - 
er 



2e' 



4L>2 



(5) 



At r > l.llD indirect magnetoexcitons attract, and at 
r < 1.1 11? they repel. The minimum of potential energy 
U{r) locates at r = « 1.671? between indirect exci- 
tons. So at large D magnetobiexciton levels correspond 
to the two-dimensional harmonic oscillator with the fre- 
quency Lu = 0.88eV(mseZ?3): £;„ = -OMe'^/{eD) + 

2V2E0 {r*/Df'^ {n + 1) where Eq = niBe"^ / {fi^e), r* ^ 
?i^e/(2m_Be'^), e'^ = 0.88e^. In the ground state the char- 
acteristic spread of magnetobiexciton aj, along QWs/GLs 
(near the mean radius of magnetobiexciton tq along 

wells/GLs) is: at = [2h/{mBUj)]^^'^ = {8r*W'^D^/^ = 
IMaeo, , where a^^ = (8r'=^)i/4D3/4. ^ex ^ ri^^l{2mBe^) 
is the two-dimensional effective Bohr radius with the ef- 
fective magnetic mass m^. Hence, the ratio of the bind- 
ing energies of the magnetobiexciton and magnetoexci- 
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FIG. 2: Dependence of temperature of the Kosterlitz- 
Thouless transition Tc — Tc{B) for the superlattice consisting 
of QWs for GaAs/AlGaAs: e — 13; and for GLs separated 
by the layer of 5*102 with e — 4.5 on the magnetoexciton den- 
sity n at D = lOnm at different magnetic fields. The solid, 
dashed and thin solid curves for QWs, dotted, dashed-dotted 
and thin dotted curves for GL at B: B = 20T, B = 15T and 
B — lOT respectively. 

ton is Ebex/Eex = 0.04 ^ 1 aX D ^ Uex, and the ratio 
of radii of the magnetoexciton and magnetobiexciton is 
ro/aex = 0.67(8r^^)i/4i:)-i/4 ^ i. So the adiabatic con- 
dition is vahd. 

The magnetobiexciton considered above has only non- 
zero quadrupole moment Q = ieD^ (the large axis of 
the quadrupole is normal to QWs/or GLs). So indirect 
magnetobiexcitons interact at long distances i? ^ Z) as 
parallel quadrupoles: U{R) = 9e'^D'^ /{eR^). 

In a 2D system, superfluidity of magnetobiexcitons ap- 
pears below the Kosterlitz-Thouless transition tempera- 
ture Tc = 7ms(r)/(2TO^) ins{T) is the density of the 



superfluid component), where only coupled vortices are 
present^. The dependence of Tc on the density of mag- 
netoexcitons at different magnetic field B for superlattice 
consisting of QWs and GLs was calculated in the ladder 
approximation for the dilute magnetobiexciton gas, and 
the results are represented on Fig. [2l 

Thus, it is shown that the low-density system of in- 
direct magnetoexcitons in a slab of superlattice of sec- 
ond type or consisting of QWs or GLs in high magnetic 
field occur to be instable due to the attraction of magne- 
toexcitons with opposite dipoles at large distances. Note 
that in spite of both QW and graphene realizations repre- 
sented by completely different Hamiltonians, the effective 
Hamiltonian in a strong magnetic field was obtained to be 
the same. Moreover, for TV excitons we have reduced the 
number of the degrees of freedom from 2 x 2 to x 2 by 
integrating over the coordinates of the relative motion of 
e and h. The instability of the ground state of the system 
of interacting two-dimensional indirect magnetoexcitons 
in a slab of superlattice with alternating electron and 
hole layers of both QWs and GLs in high magnetic field 
is claimed due to the attraction between the indirect exci- 
tons with opposite directed dipole moments. The stable 
system of indirect quasi-two-dimensional magnetobiexci- 
tons, consisting from indirect excitons with opposite di- 
rected dipole moments is considered. The low-density 
system of indirect magnetobiexcitons in superlattices is 
stable due to the quadrupole-quadrupole repulsion. So 
at the pumping increase at low temperatures the exci- 
tonic line must vanish and only magnetobiexcitonic line 
survives. 
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